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We examine the dynamics and stripe formation in a system with competing short and long range 
interactions in the presence of both an applied dc drive and quenched disorder. Without disorder, 
the system forms stripes organized in a labyrinth state. We find that, when the disorder strength 
exceeds a critical value, an applied dc drive can induce a dynamical stripe ordering transition to a 
state that is more ordered than the originating undriven, unpinned pattern. We show that signatures 
in the structure factor and transport properties correspond to this dynamical reordering transition, 
and we present the dynamic phase diagram as a function of strengths of disorder and dc drive. 
PACS numbers: 64.60.Cn, 89.75.Kd, 73.20.Qt, 71.45.Lr 



Recently, there has been considerable interest in the 
dynamics of driven elastic media in the presence of 
quenched disorder and an applied drive. Physical sys- 
tems that fall into this category include vortex lattices in 
disordered superconductors Jl|-^|, sliding charge-density 
waves H , driven electron crystals in the presence of ran- 
dom impurities (To) , sliding friction , and domain wall 
motion jl2) . One of the most studied systems in this class 
is a vortex lattice driven by an applied current through 
disordered superconductors. Higgins and Bhattacharya 
JjJ used transport measurements to map out a dynam- 
ical phase diagram for driven vortex matter based on 
transport signatures, and proposed the existence of three 
dynamical phases: a low drive pinned phase where the 
vortices do not move, a plastic phase where inhomoge- 
neous flow and tearing of the highly disordered vortex 
lattice occurs, and an elastic flow regime where the lat- 
tice slides as a whole Koshelev and Vinokur |Q in- 
vestigated the driven vortex lattice system theoretically 
and numerically, observed three phases as a function of 
increasing applied drive, and found that the disordered 
lattice in the plastic flow regime can undergo a striking 
dynamical freezing transition in which the vortex lattice 
reorders at high drives. 

Further theoretical work has shown that the recrys- 
tallized state is not fully ordered but is still strongly af- 
fected by transverse modes from the pinning. Thus the 
reordered state may form a moving smectic lattice with 
anisotropic ordering, where the vortices move in one- 
dimensional (ID) partially coupled channels aligned with 
the drive |3|-|5[] . This reordering transition to an aligned 
moving smectic state has been experimentally confirmed 
by transport || and direct imaging Q experiments. In 
addition, numerical work has confirmed the presence of 
a field-driven plastic to ordered or elastic transition ||. 
Dynamical reordering has also been studied in sliding 
charge density waves Q and driven Wigner crystals [[l0| . 

A natural question is whether these dynamical phases 
and reordering transitions can occur in other systems 
which do not form triangular lattices in the absence of 
quenched disorder. For example, many systems form 
"stripe" or "labyrinth" states including diblock- 



copolymers, magnetic domain walls pf], flux in type-I 
superconductors, water-oil mixtures JL5| , and charge or- 
dered or electron liquid crystal states in 2D electron gases 
fl6| , [r7f or superconductors |18| . These stripe- like phases 
may be disordered or destroyed by quenched disorder. 
With an applied drive, however, it may be possible to 
return to a partially ordered state. One intriguing possi- 
bility is that the reordered state may be fully aligned with 
the drive, meaning that the moving ordered state would 
be more ordered than the stationary state observed with- 
out quenched disorder. This effect could be useful for 
straightening a labyrinth forming system when ID arrays 
of aligned stripes are the desired pattern. 

Here, we model a labyrinth or stripe forming system 
in 2D by conducting molecular dynamics simulations of 
interacting particles that have a long-range Coulomb re- 
pulsion and a short range exponential attraction. Such 
a system has previously been shown to produce stripe, 
bubble and crystalline phases depending on the particle 
density and the relative strength of the attractive interac- 
tion |l8|,|l9| . We use overdamped dynamics which should 
be appropriate for systems such as colloids and magnetic 
domains. We do not take into account possible hydrody- 
namic effects. The equation of motion for a particle i is 
f = fy + f p + id = rjVi . Here, is the particle velocity, r] 
is the damping term which we set to unity, and the force 
from the other particles is = — J2j VU(r), where 

U(r) = - - Be- Rr . 
r 

For small r the repulsive Coulomb term dominates. To 
avoid a divergence of the Coulomb term we cut off the 
force at very short distances (r < 0.1) which does not 
affect any of the results presented here. We also use 
a numerically efficient summation method for the long 
range Coulomb interaction j2(J. The force from the 
quenched disorder f p comes from N p randomly placed 
attractive parabolic pins of maximum strength f p and 
radius r p = 0.3. The driving term is applied in small 
increments of 0.02 up to Fd — 6.0. We measure the re- 
sulting particle velocities V x = v i ' x by averaging the 
velocities at each fd increment value for 1000 time steps. 
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FIG. 1. (a,c,e) Real space images of the particle loca- 
tions; (b,d,f) corresponding structure factor S(k) for: (a,b) 
The moving labyrinth phase for f p = 0.5. (c,d) The intersti- 
tial flow regime for f p = 6.5. (e,f) The moving stripe regime 
for f p = 6.5. In each case the number of particles Ni = 1280 
and pinning sites N p = 400. 

We fix B = 1.25 and choose a particle density n.; = 
0.128 at which the equilibrium state in the absence of 
pinning is a labyrinth phase, rather than oriented stripes, 
which is consistent with the idea that stripe forming sys- 
tems contain self-generated randomness |ll|. The initial 
state is prepared by annealing from a high temperature 
molten state to T = 0. We conduct a series of simula- 
tions for varied f p and N p , with the number of particles 
Ni = 1280 and the system size L = 100. We have also 
run simulations for fixed n.; at various L and Ni, and find 
that only the transient responses are affected. 

We first consider the case for fewer pinning sites than 
particles, N p = 400. After annealing, we observe a dis- 
torted labyrinth pattern. Upon applying a drive, we find 
that there is no pinning threshold within our resolution 
and that the motion initially occurs by interstitial flow of 
unpinned particles around the particles that are trapped 
at the pinning sites. The overall velocity V x is still less 



than the free flow velocity since a portion of the par- 
ticles remain pinned. As the drive increases, the sys- 
tem becomes increasingly disordered and signatures of 
the labyrinth ordering are further lost. The flow in this 
regime consists of river structures where the particles flow 
in individual channels, similar to those observed in vor- 
tex systems [||. Increasing the drive, we find that the 
system can reorder into two distinct moving states: (i) 
For weak disorder, the system reorders into a labyrinth 
state similar to the pattern that forms if the system is 
annealed and without quenched disorder, (ii) For suffi- 
ciently strong disorder, the particles reorder into a stripe 
state that is aligned in the direction of the drive. In Fig. 1 
we show the real space images and the corresponding 
structure factor S(k) for these phases. 

In Fig. 1(a), the moving labyrinth phase shows stripe 
like features which have no long-range orientation. The 
structure factor has an inner ring of width correspond- 
ing to the spacing between the stripes. The ring shape 
reflects the fact that the stripes do not have a partic- 
ular orientation. The outer ring in S(k) corresponds 
to the length scale of the lattice of individual particles 
that forms within each stripe. This tendency to form a 
lattice is reflected by some remnant of six-fold order in 
the outer ring. We find the same S(k) characteristics 
for both the moving and stationary labyrinth states. In 
Fig. 1(c), we show the real space image of the particles in 
the interstitial flow region. Here the system is much more 
strongly disordered than the labyrinth regime. There are 
areas where the particles form moving channels, one par- 
ticle wide, aligned in the direction of drive. The struc- 
ture factor (Fig. 1(d)) is much more diffuse, reflecting 
the stronger disordering. There are, however, now two 
peaks for small k, which is due the smectic ordering of 
the chain structures. In Fig. 1(e) we show the high drive 
state for strong disorder, where an ordered stripe state 
forms aligned in the direction of the drive. The corre- 
sponding structure factor in Fig. 1(f) shows two peaks 
for small k corresponding to the smectic ordering of the 
stripes, in contrast to the ring structure observed in the 
labyrinth phase. Additionally the outer ring shows con- 
siderable six-fold modulations, due to the individual par- 
ticles within the stripe forming a crystal structure. 

The oriented stripe state for strong disorder is con- 
sistent with studies for vortex matter in 2D where, for 
strong disorder, the vortices form a smectic state with 
ordering occurring in only the transverse direction. For 
weaker disorder the vortex lattice anisotropy is reduced 
and the moving vortex lattice resembles the unpinned 
equilibrium state. The stripe ordering can also be viewed 
as a shear induced ordering. Shear forces are known to 
cause alignments into smectic phases pjj . In our case 
there is no global shear, but there can be local shear if a 
portion of the particles remain pinned while other parti- 
cles slide past. For weak disorder the local shear forces 
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FIG. 2. (a) The velocity V x vs F d curve (dotted line), 
and dVx/dFd curve (solid line) for f p = 6.5, N p = 400, and 
Ni — 1280. (b) The corresponding sixfold coordination frac- 
tion Pa vs Fd- Inset: the dV x /dFd curve for the same system 
as in (a), for a lower pinning strength, f p = 0.5. 

are too weak to cause a realignment of the stripes. 

We next show that the onset of the different phases 
strongly affects the transport properties. Additionally 
we show that a convenient order parameter is the frac- 
tion Pq of six-fold coordinated particles. For a perfect 
triangular lattice P§ — 1.0. We find that the stripe and 
labyrinth phases have distinct Pq values near 0.6 ± 0.015 
and 0.5±0.3, respectively, while the disordered phase has 
Pq < 0.5. In Fig. 2(a) we plot the velocity V x vs Fd for 
a sample with strong disorder, f p — 6.5, along with the 
corresponding dV x /dFd curve (solid line). In Fig. 2(b) we 
show the corresponding Pg vs Fd- There is a prominent 
peak in dV x /dFd near Fd = 1.4 and a second prominent 
peak near Fd = 2.5, each of which corresponds to a rise 
in Pq. For Fd < 1.5 the system is disordered and only 
interstitial particles move through individual channels. 
For 1.5 < Fd < 2.5, the system is still highly disordered, 
as seen from the low value of Pq « 0.35. The flow now 
occurs throughout the sample, with individual particles 
intermittently escaping briefly from pins, and then being 
rcpinncd. The first peak in dV x /dFd corresponds to the 
increase in V x caused by the previously pinned particles 
dcpinning and taking part in the motion. The additional 
small peaks in dV x /dFd correspond to specific portions 
of the pinned particles breaking free and becoming mo- 
bile in ID chains. Similar small peaks in dV x /dF curves 
observed for driven vortices are believed to correspond 
to the depinning of large clumps of vortices |Jy§]. The 
peak near Fd = 2.5 corresponds to the initial formation 
and alignment of the stripe state, which is also seen as a 
sharp rise in P§ to around 0.6. The stripe phase is fully 
formed for Fd — 2.5 where P§ plateaus. The peak in 
dV x /dF appears because the effectiveness of the pinning 
decreases as the stripes align. The dV x /dFd curve 
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FIG. 3. The dynamic phase diagram for Fd vs f p at fixed 
Ni — 1280 and N p = 400. Four phases are delineated: the in- 
terstitial flow, disordered flow, moving labyrinth, and moving 
stripe phases. 

saturates to a constant value, reflecting free flow ohmic 
behavior. We find the same features for other simula- 
tions for 2.5 < f p . For weaker pinning where a mov- 
ing labyrinth forms, the dV x /dFd curves shows only one 
peak [shown in the inset of Fig. 2(b)] corresponding to 
the initial formation of the moving labyrinth. For driven 
vortex lattices only one prominent peak is observed in 
the dV x /dF curves when the vortex lattice reorders and 
makes a transition from the plastic flow phase to the 
ohmic regime. We also find considerable hysteresis in 
the V x — Fd curves in the oriented stripe phase, while 
there is no hysteresis for the system with weak disorder. 

In Fig. 3 we show the dynamic phase diagram for Fd 
vs f p for Np — 400. We delineate four phases: the inter- 
stitial phase, determined from the dV x /dFd curves as the 
region before the first peak; the disordered flow region 
which is between the first peak in the dV x /dFd curve and 
the plateau in Pq] the moving stripe phase which is iden- 
tified from the real space images and from P$ > 0.57; and 
the moving labyrinth phase, identified from the real space 
images and from Pg w 0.5. Figure 3 shows that the ori- 
ented stripe phase appears only above a critical disorder 
strength. In Fig. 4(a) we show Pq at a fixed high drive 
of Fd = 5.0 for a series of simulations at the various f p 
used in creating Fig. 3. Fig. 4(a) clearly shows a jump 
near f p = 2.5 indicating the transition from the moving 
labyrinth to the moving stripe state, as a function of f p . 

We have also performed simulations at fixed f p and 
varied N p , and find the same general phases. The moving 
stripe phase can be still realized with as few as N p = 50 
pins. For N p < 50, the flowing state consists of coexisting 
aligned stripes and labyrinths. In Fig. 4(b) we show the 
V x — Fd curves for varied N p of 50 to 2400. The onset of 
the different phases as a function of drive changes little 
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FIG. 4. (a) The sixfold coordination fraction P& vs F p at 
a fixed value of F d = 5.0, N, = 1280 and N p = 400 show- 
ing a crossover from the moving labyrinth phase to the mov- 
ing stripe phase, (b) The velocity V x vs Fd curve for fixed 
1280 for varied N p . From top to bottom, 



F p = 6.5 and N { 

N v 



50, 400, 600, 1400, and 2400. 



with increasing N p . The value of V x in the interstitial 
region continuously decreases since more particles can be 
pinned. At very high pin density, the interstitial region 
is lost and is replaced with a pinned region. The corre- 
sponding phase diagram would have the interstitial flow 
regime replaced by a pinned flow regime. For the high 
drives F4 > 2.6 all the curves become ohmic in the free 
flow regime. We have also simulated systems with weaker 
pinning f p = 1.0 and various N p and find that the high 
drive state again becomes a moving labyrinth phase. 

In summary, we have investigated the dynamic phases 
in a driven system with quenched disorder that forms a 
disordered stripe or labyrinth phase in equilibrium with- 
out pinning. For weak pinning the system is disordered 
at low drives and reorders into a moving labyrinth phase 
at higher drives. Above a critical disorder strength, how- 
ever, the system can reorder into a moving aligned stripe 
state. This stripe state is more ordered than the unpinned 
equilibrium state. Our results confirm that the dynamical 
reordering phenomenon studied in vortex lattices can be 
applied to other systems that do not form a triangular 
lattice in equilibrium. We have shown how the phases 
can be identified through transport properties and struc- 
ture factors, and have mapped out the dynamical phase 
diagram as a function of disorder strength. Our results 
suggest that quenched disorder or other "obstacles" may 
be a useful way to engineer aligned domains in a labyrinth 



forming system. Additionally our results will be useful 
for the interpretation of transport studies in, e.g., driven 
electron-liquid crystal systems. 
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